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@) PROCEDE ET DISPOSITIF NON INVASIF DE FOCALISATION D'ONDES ACOUSTIQUES 

(§) Proc6d6 de focalisation d'ondes acoustiques utilisable 
pour obtenir une image d'un champ a observer dans un mi- 
lieu h6t6rogene dissipatif (2, 3) autour duquei on fixe des 
transducteurs acoustiques (T1-Tn, T1-Pm) formant un re- 
seau d'imagerie et un reseau cible. On suit une etape d'ap- 
prentissage au cours de laquelle on mesure des reponses 
impulsionnelles du milieu entre chaque transducteur (Ti) du 
reseau d'imagerie (5) et plusieurs transducteurs (Tj) du r6- 
seau cible (6). On en deduit des signaux de reference k 
emettre par les transducteurs du reseau d'imagerie pour 
produire une impulsion acoustique focalisee en chaque 
transducteur du reseau cible, puis de proche en proche, on 
determine des signaux de reference & emettre pour focali- 
ser une impulsion acoustique sur des points predetermines 
dans le milieu. Ces signaux de reference sont memorises et 
utilises ensuite pour generer une image acoustique du mi- 
lieu. 
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Procede et dispositif non in vasif de f ocalisation d' onde* 
acous tigues . 

La presente invention est relative aux precedes et 
dispositifs non invasif s de focalisation d'ondes 
acoustiques, notamment ultrasonores . 

Plus particulierement, 1 ' invention concerne un 
procede non invasif de focalisation d'ondes acoustiques 
dans un milieu heterogene dissipatif comprenant un milieu 
sensiblement homogene (par exemple, le cerveau) entoure au 
moins partiellement par une couche aberratrice dissipative 
(par exemple, le crane) qui genere des aberrations dans la 
propagation des ondes acoustiques, les ondes acoustiques 
etant emises depuis l'exterieur de la couche aberratrice et 
15 focalisees dans le milieu sensiblement homogene. 

Les precedes de ce type qui sont couramment 
utilises ne permettent pas d'obtenir une bonne focalisation 
des ondes acoustiques a 1'interieur du milieu, et lorsque 
ces precedes sont utilises dans des applications 
d'imagerie, ils ne permettent done pas d'obtenir une bonne 
resolution et un bon contraste d' image lorsque les 
aberrations de propagation sont importantes, par exemple 
lorsqu'on realise une echographie du cerveau depuis 
l'exterieur du crane. 

25 La P r esente invention a notamment pour but de 

pallier cet inconvenient. 

A cet effet, selon 1' invention, un proc6de de 
focalisation du genre en question est caracterise en ce 
qu'il comporte les etapes suivantes : 

30 (a) une 6ta P e initiale de positionnement au cours 

de laquelle on fixe un nombre t superieur a 2 de 
transducteurs acoustiques dans des positions predetermines 
a l'exterieur de la couche aberratrice, ces transducteurs 
etant en contact (direct ou indirect) avec ladite couche 

35 aberratrice et formant au moins : 
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- un. reseau d' imagerie qui regroupe un nombre n 
compris entre 1 et t desdits transducteurs, 

et un reseau cible qui regroupe un nombre m 
compris entre 1 et t desdits transducteurs (ces deux 
5 reseaux peuvent etre entierement distincts, ou comporter 
certains transducteurs communs, ou encore comporter chacun 
la totalite des transducteurs susmentionnes) , 

(b) une etape d' apprentissage comprenant elle-meme 
les sous etapes suivantes : 

10 (bl) une sous-etape d' apprentissage de la 

focalisation du reseau d' imagerie sur le reseau cible, 
sous-etape au cours de laquelle : 

(bll) on determine des reponses 

xmpulsionnelles hri(t) du milieu heterogene dissipatif 
15 respectivement entre chaque transducteur i du reseau 
d' imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice en correspondance respective avec des 
transducteurs du reseau cible (cette determination peut 
Stre faite par mesure directe si l'on fait emettre des 
20 impulsions acoustiques par les transducteurs du reseau 
cible, ou le cas echeant par mesure et calcul si l'on fait 
emettre des signaux acoustiques autres que des impulsions 
par les transducteurs du reseau cible, les valeurs ainsi 
mesurees et/ou calculees pouvant ensuite etre le cas 
echeant corrigees par repropagation numerique pour simuler 
des transducteurs situes directement au contact de la 
couche aberratrice si les transducteurs ne sent pas au 
contact direct de la couche aberratrice), ces reponses 
xmpulsionnelles etant memorisees sous forme numerique avec 
un certain echantillonnage temporel qui determine un nombre 
P de composantes f requentielles de l a reponse 
impulsionnelle, de frequences respectives a>k, i et ant un 
indice compris entre 1 et n qui designe un transducteur du 
reseau d' imagerie, r etant un indice compris entre 1 et m 
qui designe un point de focalisation correspondent a un 
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transducteur du reseau cible et k etant un indice compris 
entre 1 et p qui designe une composante frequent ielle, 

(bl2) a partir de ces reponses 

impulsionnelles, on calcule, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la paroi 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les differents transducteurs i du reseau d' imagerie 
genererait un signal predetermine (par exemple une 
impulsion acoustique) focalise sur le point de focalisation 
r, 

(b2) une sous-etape de focalisation en un nombre 
R de points de focalisation predetermines situes dans le 
milieu sensiblement homogene, d' indices q compris entre m+1 
et m+R, cette sous-etape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
a pas des points de focalisation lam correspondant aux 
transducteurs du reseau cible, des signaux de reference 
e'i(t,q) a faire emettre par les differents transducteurs i 
du reseau d' imagerie pour generer une impulsion focalisee 
sur ledit point de focalisation q, les signaux de reference 
e'i(t,q) etant determines pour chaque point de focalisation 
q en procedant comme suit : 

(b21) une premiere estimation de e'i(t,q), 
pour i allant de 1 a n, est calculee a partir d'au moins un 
signal de reference e'i(t,qO), qO etant 1 1 indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
ondes acoustiques dans le milieu sensiblement homogene, 

(b22) on fait emettre par les transducteurs 
du reseau d' imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e'i(t,q), 
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puis on capte avec les memes transducteurs des signaux 
Si(t,q) retrodif fuses par le milieu heterogene dissipatif 
puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la maniere suivante : 

e'iiO -► ctiM.e'iO-Tiiq)) 
ou les valeurs a i( q) et Ti (q) sent un facteur d' amplitude 
et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses , 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine, 

(b3) . on Cerise les signaux de reference 
e'i(t,q), au moins pour q compris entre m+1 et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respect ivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'^t^q), i 
etant un indice compris entre 1 et n designant 
transducteur du reseau d'imagerie. 

Grace a ces dispositions, on s'affranchit des 
aberrations de propagation des ondes acoustiques dans le 
milieu heterogene dissipatif, et on obtient une 
focalisation tres precise qui peut notamment permettre 
d'obtenir une image echographique fidele et precise d'un 
champ a observer a travers la couche aberratrice par 
retrodif fusion, lorsqu'on emet successivement des ondes 
acoustiques focalisees sur differents points du champ a 
observer et que 1'on capte les ondes acoustiques 
retrodif fusees. 

Cette focalisation precise peut egalement etre 
utilisee dans d'autres applications que 1 ' echographie, 
notamment : 

l'imagerie doppler couleur, 
- les methodes d'imagerie par elastographie, 
telles que celle decrite dans le document WO-A-00/55 616, 

les methodes d'imagerie non lineaire ("harmonic 
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imaging") , 

- les methodes de traitement par destruction 
localisee d' une partie du milieu heterogene dissipatif, 
notamment par hyperthermie, 

les methodes de mesure de parametres 
d' absorption optique des tissus avec activation par 
ultrasons, etc. 

Dans des modes de realisation preferes de 
1' invention, on peut eventuellement avoir recours en outre 
a l'une et/ou a 1' autre des dispositions suivantes : 

- au cours de la sous-etape (bll) , lorsqu'au moins 
certains transducteurs (du reseau cible et/ou du reseau 
d'imagerie) sont en contact avec un milieu homogene 
intermediaire (par exemple un gel) lui-meme en contact avec 
la couche aberratrice, on corrige les reponses 
impulsionnelles hri(t) par repropagation numerique pour 
simuler des transducteurs situes directement au contact de 
la couche aberratrice ; 

- la sous-etape (bl2) comporte elle-meme les sous- 
etapes suivantes : 

(bl21) on determine p matrices de transfert 
H(cok) = [Hri(a>k)], i allant de 1 a n et r allant de 1 a m, 
ou Hri(cok) est la valeur, a la frequence cok, de la 
transformee de Fourier de la reponse impulsionnelle hri(t), 

(bl22) on determine pour chaque point de 
focalisation r correspondent a un transducteur du reseau 
cible, n composantes Ei(e>k,r), i variant entre 1 et n, 
telles que F (cok, r ) =H (cok) .E(cok,r) , ou E(tok,r) = [Ei(tok,r)] 
est un vecteur a n composantes, F(cok,r) est un vecteur a m 
composantes Fl(cok,r), 1 variant entre 1 et m, ces m 
composantes Fl(cok,r) correspondant a une focalisation 
souhaitee des ondes acoustiques a la frequence cok sur le 
point de focalisation r correspondant a un transducteur du 
reseau cible, 
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(bl23) on en deduit, pour chaque point de 
localisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels e (t, r) = [ei (t , r) ] , 

i variant entre 1 et n, ou «(r,r)=^£/(^,r).^' en notation 

complexe, ces signaux ei(t,r) etant adaptes pour que 
1' emission de ces derniers respectivement par les 
differents transducteurs i du reseau d' imagerie genere une 
impulsion acoustique focalisee sur le point de focalisation 
r du reseau cible, 

(bl24) une sous-etape de correction des 

aberrations generees par la couche aberratrice entre le 
milieu sensiblement homogene et chaque transducteur cible 
r, ces aberrations etant estimees sur la base des mesures 
precedemment effectuees, les aberrations ainsi estimees 
etant utilisees pour calculer lesdits signaux temporels de 
reference e'i(t,r) ,- 

- au cours de la sous-etape (bl22) on calcule p 
matrices H" 1 (cok) , respectivement par regularisation et 
inversion des matrices de transfert H(cok), et pour chaque 
transducteur r du reseau cible, on calcule le vecteur 
E(Q)k,r) par la formule : 

Efo^rJ^IT^oak) .F(cak, j) ; 

- au cours de 1 • etape (bl22) , les' composantes 
Fl((ok,r) du vecteur F(<Dk,r) correspondant a la repartition 
spatiale du champ desire a la frequence cok, sont egales a 
0 pour l*r et egale a 1 pour l=r ; 

- au cours de la sous-etape (bl24), on assimile la 
parol aberratrice, au voisinage de chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, a un filtre a reponse impulsionnelle finie, defini a 
chaque frequence cok par une amplitude Gr(cok) et une phase 
<Mcok), la sous-etape (bl24) comportant elle-meme les sous- 
etapes suivantes : 
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(bl241) on calcule, pour chaque frequence cok, 
1' amplitude Gr (cok) et la phase <(> r (cok), a partir soit des 
signaux ei(t,r), soit des vecteurs E (cok, r) , 

(bl242) on calcule p matrices de transfert 
corrigees H' (cok) = [H'ji(cok)], ou 

(bl243) on determine pour chaque transducteur r 
du reseau cible, n composantes E'i(cok,r), i variant entre 
1 et n, telles que F (cok, r ) =H / (cok) . E' (cok, r ) , ou E' (cok, r ) = 
[Ei(tok,r)] est un vecteur a n composantes, F(cok,r) est un 
vecteur a m composantes Fl(cok,r), 1 variant entre 1 et m, 
ces m composantes Fl(cok,r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 
cok sur le point de focalisation r correspondant a un 
transducteur du reseau cible, 

(bl244) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels de reference 
e' (t , r) = [e' i (t, r) ] , i variant entre 1 et n, ou 

e\(t,r)=^Ei((Dk 9 r)^? a * t en notation complexe ; 

- au cours de la sous-etape (bl241) , on calcule 
1' amplitude Gr(cok) et la phase <|> r (cok) comme suit : 



Gr(<Dky 



1 



Ei(cok 9 r)Ei (6>*,r) 



ou 



Ei* est la valeur complexe conjuguee de Ei, 

et Ax(rO # r, i)=(d(rO, i) -d(r, i) ) /c, d(r,i) etant la 
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distance entre le transducteur i et le point de 
focalisation r, et d(r0,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO ; 

- la sous-etape (bl2) comporte elle-rneme les sous- 
etapes suivantes : 

(bl21) on determine p matrices de transfert 
H(a>k) = [Hri(a>k)], i allant de 1 a n et r allant de 1 a m, 
ou Hri(cok) est la valeur, a la frequence cok, de la 
transform^ de Fourier de la reponse impulsionnelle hri(t), 

(bl22') on corrige les matrices de transfert 
H (a>k) pour s'affranchir des aberrations generees par la 
paroi aberratrice au voisinage de chaque point de 
focalisation r, cette correction etant effectuee a partir 
des reponses impulsionnelles hri(t) precedemment 
15 determines, et on obtient ainsi des matrices de transfert 
corrigees H' (cok) , 

(bl23' ) on determine pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, n composantes E'i(cok,r), i variant entre 1 et n, 
20 telles que F (cok, r ) =H' (cok).E' (cok, r ) , ou E'«ok,r) = 
[E'i(ok,r>] est un vecteur a n composantes, F((ok,r) est un 
vecteur a m composantes Fl (cok, r ) , 1 variant entre 1 et m, 
ces m composantes Fl (cok, r) correspondant a une 
focalisation souhaitee des ondes acoustiques a la frequence 
25 cok sur le point de focalisation r correspondant a un 
transducteur du reseau cible, 

(bl24') on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
C J- ble ' Un vecteur de n signaux temporels 

e' (t,r) = [e'i(t,r) ] , i variant entre 1 et n, ou 

*(ts^Ei{G>k,r)e J4 *' en notation complexe, les signaux 
e'i(t,r) etant lesdits signaux de reference ; 

au cours de la sous-etape (bl23' ) on calcule p 
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matrices H'-(a>k), respect ivement par regularisat ion et 
inversion des matrices de transfert H' (wk) , et pour chaque 
transducteur r du reseau cible, on calcule le vecteur 
E'(cok,r) par la formule : 

E' (cok, r) =H' - 1 (cok) . F (ok, j ) ; 

- au cours de 1 • etape (bl23' ) , les composantes 
Fl(cok,r) du vecteur F(cok,r) correspondant a la repartition 
spatiale du champ desire a la frequence cok, sont egales a 
0 pour l*r et egale a 1 pour l=r ; 

- au cours de la sous-etape (bl22' ) , on assimile 
la parol aberratrice, au voisinage de chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, a un filtre a reponse impulsionnelle finie, defini a 
chaque frequence cok par une amplitude Gr (cok) et une phase 
•Mcok), la sous-etape (b!22' ) comportant elle-meme les 
sous-6tapes suivantes : 

(bl22'l) on calcule, pour chaque frequence 
cok, 1' amplitude Gr(cok) et la phase +r «ok), a partir des 
reponses impulsionnelles precedemment determines, 

(bl22'2) on calcule p matrices de transfert 
corrigees H' (cok)=[H'ji(cok) ] , ou 

G A*>k) 

- au cours de la sous-etape (bl22'l) on calcule, 
pour chaque frequence cok, 1' amplitude Gr (cok) et la phase 
25 <(» r (cok) / de la mani^re suivante : 



J]?Hri(G> k )JI rt \a> k ) 



ou 
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H*ri designe la valeur complexe conjuguee de Hri, 
et AT(rO,r,i) = (d(rO / i)-d(r,i))/ C/ d(r,i) etant ' la 
distance entre le transducteur i et l e point de 
focalisation r, et d(r0,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO ; 

- au cours de l'etape (c) , on fait suivre la sous- 
etape (cl) par les sous-etapes suivantes : 

(c2) on fait capter par lesdits transducteurs du 
reseau d' imagerie des signaux s i( t) retrodif fuses par le 
milieu heterogene dissipatif, 

<c3) on convolue le signal de reference emis par 
chaque transducteur du reseau d' imagerie avec le signal 
retrodiffuse capte par ce transducteur, 

(c4) puis on somme les produits de convolution 
ainsi obtenus, 

l'etape (c) etant renouvelee pour une pluralite de points 
situes dans le milieu sensiblement homogene ; 

- au cours de la sous-etape (b21) , l a premiere 
estimation de chaque signal de reference est e'i(t,q) = 
e'i(ts+ei(q),q 0 ) pour chaque point de focalisation q, q 0 
etant 1 • indice d'un point de focalisation proche du point q 
pour lequel le signal de reference a deja ete determine, 
ei(q) etant un retard egal a une valeur 5i (q) / c , ou c est 
la celerite moyenne des ondes acoustiques dans le milieu, 
et 6i(q) est egal a une difference entre d' une part, une 
distance entre le transducteur i du reseau d' imagerie et le 
point de focalisation q0, et d' autre part, une distance 
entre le transducteur i du reseau d' imagerie et le point de 
focalisation q, 

- au cours de la sous-etape (b2), lorsqu'au moins 
certains transducteurs d' indice v du reseau d'imagerie ne 
sont pas directement au contact de la couche aberratrice, 
on corrige les signaux e' v (t,q) correspondants par 
repropagation numerique pour simuler des transducteurs 
places en contact direct avec la couche aberratrice ; 
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- au cours de la sous-etape (b22), on recherche 
les valeurs a^q) et ti(q) pour maximiser le critere de 
coherence C suivant : 



C = 





n 


2 


< 


M 


> 



ou 



"-E<ks,('-*7,?)| 2 > 



8i(*><l) = s i (t)<g)e i (t,q) , ® representant 1' operation de 
convolution, 

et <> represente une moyenne temporelle ; 

- au cours de la sous-etape (b22) , les valeurs 
Ti(q) sont calculees en maximisant une fonction 
d'intercorrelation, pour des transducteurs voisins du 
reseau d'imagerie, des signaux g L (t , q) et g i+1 (t,q) ; 

- au cours de la sous-etape (b22) , les valeurs 
cti(q) sont calculees de maniere a egaliser sur l'indice i 
1' amplitude maximale des fonctions gi(t,q) ; 

- au cours de la sous-etape (b22) , les valeurs 
ai(q) et x^q) sont calculees en realisant une 
intercorrelation, pour des transducteurs voisins du reseau 
d'imagerie, des signaux g ± (t,q) et g i+1 (t,q) ; 

- au cours de la sous-etape (b22) , les valeurs 
20 ai (q) et Ti(q) sont calculees de maniere a egaliser sur 

l'indice i 1' amplitude maximale des fonctions gi(t,q) ; 

- la sous-etape (b22) relative a chaque point de 
focalisation q est realisee immediatement apres la sous- 
etape (b21) relative au meme point de focalisation q ; 

' ib ~ le m ilieu heterogene dissipatif est constitue 

par le cerveau entoure par le crane ; 

- le reseau d'imagerie et le reseau cible sont 
deux reseaux distincts disposes de part et d'autre du 
milieu heterogene dissipatif ; 

- tous les transducteurs appartiennent a la fois 
au reseau d'imagerie et au reseau cible ; 



30 



pwcnnnn. ,co ooi K 717ai i 



2815717 



15 



20 



12 



les ondes acoustiques sont des ondes 
ultrasonores. 

Par ailleurs, 1' invention a egalement pour objet un 
dispositif adapte pour mettre en oeuvre le precede defini 
5 ci-dessus. 

D'autres caracteristiques et avantages de 
1' invention apparaitront au cours de la description 
suivante d'une de ses formes de realisation, donnee a titre 
d'exemple non limitatif, en regard du dessin joint. 
10 Sur le clessin, la figure unique represente un 

dispositif d'imagerie ultrasonore selon une forme de 
realisation de 1» invention. 

Le dispositif 1 d'imagerie ultrasonore represente 
sur le dessin est adapte pour realiser une image 
echographique par ultrasons du cerveau 2 d'un patient (a 
des frequences par exemple de l'ordre de 1 a 3 MHz) , depuis 
l'exterieur du crane 3, le cerveau 2 constituant un milieu 
sensiblement homogene pour la propagation des ondes 
acoustiques et le crane 3 constituant une couche 
aberratrice dissipative, de sorte que 1' ensemble de la 
boite cranienne 2,3 constitue un milieu heterogene 
dissipatif. 

En variante, 1' invention serait applicable 
notamment : 

25 ' * l'imagerie acoustique de tout autre milieu 

heterogene dissipatif non homogene comprenant un milieu 
sensiblement homogene entoure par une couche dissipative 
relativement peu epaisse generant des aberrations dans la 
propagation des ondes ultrasonores, 

30 - ou a tout autre precede impliquant au moins une 

focalisation a 1' emission dans un tel milieu. 

Dans 1' exemple represente sur le dessin, le 
dispositif d'imagerie 1 comporte un micro-ordinateur 4, ou 
tout autre dispositif de commande et/ou de visualisation 
35 des images ultrasonores, ce micro-ordinateur comportant 
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classiquement un clavier 4a associe eventuellement a 
d'autres interfaces de commande et un ecran 4b permettant 
de visualiser les images du cerveau 2. 

Par ailleurs, le dispositif d'imagerie 1 comporte 
deux reseaux 5, 6 de transducteurs ultrasons Tl, T2... Tn et 
T'l, T'2... T'm formant par exemple deux barrettes lineaires 
de transducteurs qui sont disposees de part et d' autre du 
crane 3 de 1 ' utilisateur , dans des positions geometriques 
predetermines l'une par rapport a 1' autre, chaque barrette 
de transducteur 5, 6 etant mise en contact avec le crane 3 
par 1 • intermediate d'une couche 7 de gel ou similaire. 

Les differents transducteurs Tl, T2... Tn et T'l, 
T'2... T'm peuvent etre commandes directement par le micro- 
ordinateur 4 , ou de preference par une unite centrale 
electronique CPU contenue par exemple dans une baie 
electronique 8 et elle-meme commandee par le micro- 
ordinateur 4 . 

Avantageusement, chacun des transducteurs Tl, T2... 
Tn, T'l, T'2, T'm est relie a un echantillonneur, 
20 respectivement El, E2... En, E'l, E'2, E'm, et chaque 
echantillonneur est lui-meme relie a une memoire, 
respectivement Ml, M2... Mm, M'l, M'2... M'm et a une unite 
centrale CI, C2, ...Cm, C 1, C 2, ... Cm. Ces memoires et ces 
unites centrales sont a leur tour relies, directement ou 
indirectement, a 1' unite centrale CPU susmentionnee, qui 
est par ailleurs reliee au moins a une memoire centrale M. 

Le dispositif qui vient d'etre decrit fonctionne 
comme suit . 

Initialement, les deux reseaux de transducteurs 5, 
6 sont fixes de part et d' autre du crane 3 du patient, dans 
^lesdites positions predetermines. A cet effet, les reseaux 
de transducteurs 5, 6, dits respectivement reseau 
d'imagerie et reseau cible, peuvent etre portes par un 
support rigide tel qu'un casque (non represents dispose 
35 autour de la tete du patient. 
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Ensuite, l e dispositif suit une etape 
d'apprentissage de quelques minutes (avan.tageusement 1 a 3 
minutes) permettant de prendre en compte 1 • ensemble des 
aberrations de propagation dues au caractere non homogene 
5 du milieu dissipatif forme par le crane 3 et le cerveau 2. 

Au cours de cette etape d'apprentissage on fait 
d'abord emettre successivement par chacun des transducteurs 
Tl, T2... Ti,... Tn du reseau d'imagerie 5, une impulsion 
acoustique, et pour chaque impulsion emise par 1 • un des 

10 transducteurs Ti du r^sMn H'im^ Q ,i« 

uu reseau a imagerie, on enregistre le 

signal capte par les transducteurs T'l, T'l... T'r,... T'm du 
reseau cible 6, c'est-a-dire la reponse impulsionnelle 
hn(t) du milieu heterogene dissipatif entre le 
transducteur i considere du reseau d'imagerie 5 et chaque 
15 transducteur j du reseau cible 6. 

Chaque reponse impulsionnelle hri (t) est 
enregistree sous forme numerique avec un certain 
echantillonnage temporel qui determine un certain nombre p 
de composantes f requentielles monochromatiques de la 
20 reponse impulsionnelle, correspondant chacune a une 
frequence ok, k etant un indice compris entre 1 et p. 

Dans le cas envisage ici, ou au moins certains 
transducteurs du reseau cible et/ou du reseau d'imagerie ne 
sont pas directement au contact avec la couche aberratrice 
3, on corrige les reponses impulsionnelles pour simuler des 
transducteurs virtuels disposes au contact de ladite couche 
aberratrice. La position de la couche par rapport aux 
transducteurs pent eventuellement etre obtenue par imagerie 
conventionnelle (echographie ultrasonore, scanner X, IRM, 
30 etc.). Les reponses impulsionnelles corrigees sont 
calculees par un algorithme de repropagat ion numerique 
connu, decrit notamment dans les articles suivants : 

"Ultrasonic beam steering through inhomogeneous 
layers with a time reversal mirror", C.DORME, M. FINK, IEEE 
35 Transactions Ultrasonics, Ferroelectric and Fre q uency 
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Control, 43 (1) , janvier 1996, p 167-175, 

"Focusing and steering through absorbing and 
aberrating layers : Application to ultrasonic propagation 
through the skull" Journal of Acoustical . Society of 
America, 103 (5), Mai 1998, p. 2403-2410, 

et "Propagation and backpropagation for 
ultrasonic wavefront design" Liu, D.-L., and Waag, R. C. 
IEEE Trans, on Ultras. Ferro. and Freq. Contr. 44(1): 1-13 
(1997) . 

Dans ce qui suit, hri(t) denommera done les reponses 
impulsionnelles pour des elements (reels ou virtuels) 
situes contre la couche aberratrice. De plus, les elements 
virtuels ou r6els situes contre la couche aberratrice 3 
seront appeles ci-apres "points de f ocalisation" d' indice r 
compris entre 1 et m. 

Lorsque les transducteurs du reseau d'imagerie 5 
emettent des signaux acoustiques e ± (t), ces signaux 
g6nerent au niveau des transducteurs r du reseau cible 6 
des signaux acoustiques fr(t) s'exprimant comme suit : 

X0=Zhri(t)®ei(t), 

oil ® represente l'operateur de convolution temporelle. 

Apres transformee de Fourier, cette equation 
devient : 

F(cok) = H(ok) .E(Q)k) , ou : 

H(cok) est la matrice de transfert, de taille 
m*n, entre les transducteurs Ti du reseau d' imagerie et les 
transducteurs Tr du reseau cible : les composantes Hri (cok) 
de cette matrice sont les composantes des transformees de 
Fourier des reponses impulsionnelles hri(t) a la frequence 
<ok, 

E(ook) est un vecteur dont les composantes 
Ei(cok) sont les composantes de la transformee de Fourier 
des signaux ei(t) susmentionnes a la frequence cok, 
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et F(a>k) est un vecteur dont les composantes 
Fjttok) sont les composantes de la transformee de Fourrier 
des signaux f j (t) susmentionnes a la frequence ok. 

Par inversion de chaque matrice de transfert H (ok) , 
on peut done determiner le vecteur E(ok,j) qui est adapte 
pour generer au niveau du point de localisation r 
correspondant au transducteur T' r du reseau cible, un 
vecteur F(a>k,j) dont toutes les composantes se rapprochent 
au mieux de l'objectif initialement fixe (de preference 
toutes egales a zero, sauf la composante d'indice j 
correspondant au transducteur T • j , qui est egale a 1 
lorsqu'on veut emettre une impulsion acoustique au niveau 
du point de focalisation r) , grace a la relation :.. 

E(a>k, j)= H^tok) .F(iok, j) , 
oil H -1 (o>k) est la matrice inverse de H(cok). 

H " 1(fi)k) P eut etre calculee par exemple par 
decomposition en valeurs singulieres, ceci permettant une 
regularisation de 1' inversion de la matrice H(cok). 

Puis on determine, par transformee de Fourier 
inverse des differentes composantes Ei(ok,j) du vecteur 
E(G)k,j), les differents signaux de reference ei(t,j) qui, 
lorsqu'ils sont emis par les differents transducteurs Ti du 
reseau d'imagerie 5, sont adaptes pour focaliser une 
impulsion acoustique (ou le cas echeant un autre signal 
acoustique) au niveau du point de focalisation r. On a done 
realise une focalisation du reseau d' imagerie 5 sur chaque 
transducteur du reseau cible 6, par filtre inverse spatio- 
temporel . 

L' unite centrale CPU suit ensuite un processus 
d'apprentissage des aberrations dues a la paroi du crane 3 
au niveau du reseau cible. 

Au cours de ce processus, ces aberrations sont 
considerees comme un filtre a reponse impulsionnelle finie. 

Dans le domaine de Fourier, ce filtre est defini a 
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chaque frequence ©k une amplitude Gr(cok) et une phase 
<|>r(ook) . 

Pour calculer ces coefficients, on compare la phase 
et 1' amplitude de 1' ensemble des vecteurs Er. A cet effet, 
on commence par eliminer les dephasages introduits par les 
differences de marche entre les transducteurs d' imagerie Ti 
et les differents points de focalisation indices r. Ceci 
revient a choisir un point de focalisation particulier rO, 
et introduire pour les autres un dephasage lineaire avec la 
pulsation : exp (-jAx (rO, r, i) ©) avec At (rO, r, i) = (d (rO, i) - 
d(r,i))/c ou d(r,i) est la distance entre le transducteur i 
et le point de focalisation r et c est la vitesse moyenne 
des ondes acoustiques dans le milieu a imager, en 
1' occurrence le cerveau 2. 

Cette correction effectuee, les differences 
d' amplitude et de phase entre les vecteurs Er sont 
attributes a la couche aberratrice 3 situees contre le 
reseau cible. On calcule alors pour chaque point de 
focalisation r, le facteur de gain Gj (©k) et le facteur de 
phase <|>j(a)k) : 



f££i(^t,rO)J?,*(fi»,rt)) 

1=l£ n 



Gr(a>k)= 

'(cok,r) 



ou Ei* est la valeur complexe conjuguee de Ei . 

Ces couples {Gj (©k) , ^ (cok) } correspondent au 
facteur d' attenuation et au dephasage relatifs introduits a 
chaque frequence par la portion de couche aberratrice 3 
situee contre le point de focalisation r. lis caracterisent 
done finalement les aberrations introduites par la portion 
de couche aberratrice situee contre le reseau cible. 

On elimine ensuite les aberrations introduites par 
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la couche aberratrice 3 situee contre les transducteurs 
cibles, dans 1' ensemble des p matrices H (cok) = [H ji (cok) ] , 
def inies ci-dessus . 

Pour cela, on calcule un nouvel ensemble de 
matrices de transfert H' (cok) = [H' j i (cok) ] caracterisant la 
propagation entre le r6seau d' imagerie et le reseau cible 
dans un milieu virtuel pour lequel seules subsistent les 
aberrations situees contre le reseau d' imagerie : 

G j\<*>k) 

On determine alors pour chaque transducteur r du 
reseau cible, n composantes E'i(cok,r), i variant entre 1 
et n, telles que F (cok, r ) =H' (cok) . E' (cok, r ) , oil E' (cok, r) = 
[E'i(cok,r)] est un vecteur a n composantes, F(©k,r) est un 
vecteur a m composantes Fl(cok,r), 1 variant entre 1 et m, 
ces m composantes Fl(o>k,r) correspondant £ une 
focalisation souhait6e des ondes acoustiques a la frequence 
cok sur le point de focalisation r correspondant a un 
transducteur du reseau cible. 

On en deduit, pour chaque point de focalisation r 
correspondant a un transducteur du reseau cible, un vecteur 
de n signaux temporels de reference e' ( t , r ) = [e' i (t, r) ] , i 

variant entre 1 et n, ou ^(^r)=^/(^,r). e y<ufc/ en notation 

*-i 

complexe . 

Ces signaux de reference e'i(t,j) sont adaptes pour 
que 1' emission de ces derniers respectivement par les 
diff6rents transducteurs i du reseau d' imagerie g6ndre une 
impulsion acoustique focalisee sur le transducteur j du 
reseau cible en 1' absence de la couche aberratrice situee 
contre le reseau cible. 

On notera qu'en variante, les signaux de reference 
pourraient etre determines de la fa?on suivante apres la 
determination des reponses impulsionnelles hri(t) et des p 



2815717 



19 

matrices de transfert H(cok) : 

on corrige les matrices de transfert H(cok) pour 
s'affranchir des aberrations generees par la paroi 
aberratrice 3 au voisinage de chaque point de focalisation 
r, cette correction etant effectuee a partir des reponses 
impulsionnelles hri(t) precedemment determinees, et on 
obtient ainsi des matrices de transfert corrigees H' (cok) , 



pour chaque point de focalisation r correspondant a un 
transducteur du reseau cible, n composantes E'i(cok,r), i 
variant entre 1 et n, telles que F (cok, r ) =H' (cok) . E' (cok, r) , 
ou E'(cok,r) = [E'i(cok,r)] est un vecteur a n composantes, 
F(cok,r) est un vecteur a m composantes Fl (cok, r) , 1 variant 
entre 1 et m, ces m composantes Fl(cok,r) correspondant a 
une focalisation souhaitee des ondes acoustiques a la 
frequence cok sur le point de focalisation r correspondant 
a un transducteur du reseau cible, 



focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels e'(t,r) = 



notation complexe, les signaux e'i(t,r) etant lesdits 
signaux de reference. 

Avantageusement, lors du calcul des matrices 
H'(a>k), on assimile la paroi aberratrice, au voisinage de 
chaque point de focalisation r correspondant a un 
transducteur du reseau cible, a un filtre a reponse 
impulsionnelle finie, defini a chaque frequence cok par une 
amplitude Gr(cok) et une phase <|> r (cok) calcules comme suit : 



on determine par inversion des matrices H' (cok) , 



et on en deduit, pour chaque point de 




en 
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1=1 

H*ri designe la valeur complexe conjuguee de Hri, 
et Ax(rO,r,i) = (d(rO,i)-d(r,i) ) /c, d(r,i) etant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO. 

On calcule ensuite p matrices de transfert 
corrigees H' (cok) = [H' ji (©k) ] , ou 

Hji&k) = Hj^). g ^ e -M">y f qui S ervent a determiner les 

vecteurs E' i (cok) comme explicite precedemment et done les 
differents signaux de reference e'i(t,r), r allant de 1 a 



m. 



un 



L' unite centrale CPU procede ensuite a 
apprentissage de focalisation en un nombre R de points de 
focalisation predetermines situes dans le cerveau 2, 
d' indices q compris entre m+1 et m+R, cette sous-etape 
consistant a determiner pour chacun de ces points de 
focalisation q, en s'eloignant pas a pas des transducteurs 
du reseau cible, des signaux de reference e'i(t,q) a faire 
emettre par les differents transducteurs du reseau 
d'imagerie pour generer une impulsion focalisee sur ledit 
point de focalisation q. 

Les signaux de reference e'i(t,q) sont initialement 
determines, pour chaque nouveau point de focalisation q, 
sous la forme e'i(t,q) = e' i (t+0i (q) , q0) pour chaque point 
de focalisation q, q0 etant 1 • indice d ' un point de 
focalisation proche du point q pour lequel le signal de 
reference a deja ete determine, le retard 9i (q) etant 
initialement egal a une valeur 8i (q) /c, ou c est la valeur 
moyenne de celerite des ondes acoustiques dans le milieu, 
et 5i(q) est egal a une difference entre d' une part, une 
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distance entre le transducteur i du reseau d' imagerie et le 
point de focalisation qO, et d' autre part, une distance 
entre le transducteur i du reseau d' imagerie et le point de 
focalisation q. 

Dans le cas ou certains transducteurs d' indice v du 
reseau d' imagerie ne sont pas situes contre la couche 
aberratrice, il est en outre souhaitable de corriger les 
signaux de reference les signaux e' v (t,q) correspondants 
par repropagation numerique des transducteurs virtuels 
(situes contre la couche aberratrice 3) aux transducteurs 
reels (separes de ladite couche 3 par du gel 7 ou autre), 
de facon connue en soi, par le procede inverse de celui 
decrit ci-dessus en ce qui concerne les reponses 
impulsionnelles . 

On fait ensuite emettre par les transducteurs du 
reseau d' imagerie, par iterations, les estimations obtenues 
des signaux de reference e'i(t,q), puis on capte avec les 
memes transducteurs les signaux, Si (t,q), retrodif fuses par 
le milieu heterogene dissipatif. 

Puis, on modifie pour 1' iteration suivante ces 
signaux de reference e'i(t,q) de la maniere suivante : 

«i(0 -» <*/(?). 4 C-r,(f)) 
ou les valeurs Ti(q) et <Xi(q) sont un retard et un facteur 
d' amplitude correct if s, calculees pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses . 

Avantageusement, l e critere C pourra etre le 
suivant : 



< 

C = — 



2>/ Si ('-*-,,?) 

i=l 



OU 



n -'E < \ a i-S i (.t-r i ,qf > 

gi(t>q) = s i (t)®e i (t,q) , ® representant l'operation de 
convolution, 

et <> represente une moyenne temporelle. 
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Dans ce processus d' optimisation, les valeurs Ti(q) 
peuvent etre calculees a chaque iteration de facon a 
maximiser une fonction d' intercorrelation, pour des 
transducteurs voisins du reseau d'imagerie, des signaux 
gi(t,q) et g i+ i(t,q) susmentionnes, et les valeurs ai(q) 
peuvent etre calculees de. maniere a egaliser sur 1' indice i 
1' amplitude maximale des fonctions gi(t,q). 

On optimise ainsi les signaux de reference 
e'i(t,q), i allant de 1 a n, de facon qu'ils produisent un 
signal acoustique focalise de facon precise sur le point de 
focalisation q situe dans le cerveau. Ce processus 
d' optimisation a deja ete exj>licite plus en detail par 
Mallart et al . (The Van Cittert-Zernike theorem in pulse 
echo measurements, J. Acoust . Soc. Am. 90(5), novembre 
1991, p. 2716-2727 ,- Adaptative focusing in scattering 
media through sound speed inhomogeneities : the Van Cittert 
Zernike approach and focusing criterion, J. Acoust. Soc. 
Am. 96(6), decembre 1994, p. 3721-3732). 

Lorsque cette optimisation est terminee pour un 
point de focalisation q, par exemple apres 2 ou 3 
iterations lorsque le critere C a atteint une valeur 
predetermine (notamment proche de 2/3), on passe au point 
de focalisation q+1 suivant, etc. 

Les signaux de reference e'i(t,q) ainsi obtenus 
sont memorises, par exemple dans les memoires Ml-Mn. 

Une fois l'etape d' apprentissage terminee, on peut 
notamment realiser des images echographiques du cerveau 2, 
eventuellement a un rythme eleve pouvant atteindre la 
cadence d'un echographe standard, par exemple 20 a 30 
images par seconde. Pour realiser chacune de ces images, on 
precede comme suit, pour chaque point de focalisation q 
appartenant au champ a observer : 

on fait emettre respect ivement par les 
transducteurs Ti du reseau d'imagerie, lesdits signaux de 
reference ei(t,q), 
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- puis on fait capter par lesdits transducteurs du 
reseau d'imagerie des signaux si(t) retrodif fuses par le 
milieu viscoelastique, 

- on convolue le signal de reference ei(t,q) emis 
par chaque transducteur du reseau d'imagerie avec le signal 
retrodiffuse si(t) capte par ce transducteur, 

- puis on somme les produits de convolution ainsi 

obtenus . 

On notera que les differentes operations 
susmentionnees realisees au cours de l'etape 
d'apprentissage ou de l'etape d'imagerie peuvent etre soit 
programmees dans 1' unite centrale CPU, soit realisees en 
tout ou partie par des circuits specialises. 

Par ailleurs, on notera egalement que tous les 
transducteurs Ti, T'r pourraient etre utilises pour 
realiser les images echographiques du cerveau. Dans ce cas, 
le reseau d'imagerie serait confondu avec le reseau cible 
et chacun de ces deux reseaux comprendrait tous les 
transducteurs, le f onctionnement decrit ci-dessus 
s'appliquant alors mutatis mutandis. 
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REVINDICATIONS 

1. Procede non invasif de focalisation d' ondes 
acoustiques dans un milieu heterogene dissipatif (2,3) 
coraprenant un milieu sensiblement homogene (2) entoure au 
moins partiellement par une couche aberratrice dissipative 
(3) qui genere des aberrations dans la propagation des 
ondes acoustiques, les ondes acoustiques etant emises 
depuis l'exterieur de la couche aberratrice (3) et 
focalisees dans le milieu sensiblement homogene (2), 
caracterise en ce qu'il comporte les etapes suivantes : 

(a) une ©tape initiale de positionnement au cours 
de laquelle on fixe un nombre t superieur a 2 de. 
transducteurs acoustiques (Tl-Tn, T'l-T'm) dans des 

15 positions predetermines a l'exterieur de la couche 
aberratrice (3), ces transducteurs etant en contact avec 
ladite couche aberratrice et formant au moins : 

- un reseau d' imagerie (Tl-Tn) qui regroupe un 
nombre n compris entre 1 et t desdits transducteurs, 

20 ~ et un reseau cible (T'l-T'm) qui regroupe un 

nombre m compris entre 1 et t desdits transducteurs, 

(b) une etape d' apprentissage coraprenant elle-meme 
les sous etapes suivantes : 

(bl) une sous-etape d' apprentissage de 
25 focalisation du reseau d' imagerie sur le reseau cible, 
sous-etape au cours de laquelle : 

(bll) on determine des reponses 

impulsionnelles hri (t) du milieu heterogene dissipatif, 
respectivement entre chaque transducteur i du reseau 
d' imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice (3) en correspondance respective avec 
des transducteurs du reseau cible, ces reponses 
impulsionnelles etant memorisees sous forme numerique avec 
un certain echantillonnage temporel qui determine un nombre 
35 p de composantes f requentielles de la reponse 
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impulsionnelle, de frequences respectives cok, i etant un 
indice compris entre 1 et n qui designe un transducteur du 
reseau d'imagerie, r etant un indice compris entre 1 et m 
qui designe un point de focalisation correspondant a un 
transducteur du reseau cible et k etant un indice compris 
entre 1 et p qui designe une composante f requentielle, 

(bl2) a partir de ces reponses 

impulsionnelles, on calcule, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la parol 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les differents transducteurs i du reseau d'imagerie 
genererait un signal predetermine focalise sur le point de 
focalisation r, 

(b2) une sous-etape de focalisation en un nombre 
R de points de focalisation predetermines situes dans le 
milieu sensiblement homogene, d' indices q compris entre m+1 
et m+R, cette sous-etape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
a pas des points de focalisation 1 a m correspondant aux 
transducteurs du reseau cible, des signaux de reference 
e'i(t,q) a faire emettre par les differents transducteurs i 
du reseau d'imagerie pour generer une impulsion focalisee 
sur ledit point de focalisation q, les signaux de reference 
e'i(t,q) etant determines pour chaque point de focalisation 
q en procedant conune suit : 

(b21) une premiere estimation de e'i(t,q), 
pour i allant de 1 a n, est calculee a partir d'au moins un 
signal de reference e'i(t,q0), q0 etant 1 ■ indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
ondes acoustiques dans le milieu sensiblement homogene (2), 
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(b22) on fait emettre par les transducteurs 
du reseau d'imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e'i(t,q), 
puis on capte avec les memes transducteurs des signaux 
5 si(t r q) retrodif fuses par le milieu heterogene dissipatif, 
puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la maniere suivante : 

e)(t) aiiqyeiit-Ttiq)) 
ou les valeurs cti(q) et xi(q) sont un facteur d' amplitude 

10 et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits " signaux retrodif fuses, 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine, 

(b3) on memorise les signaux de reference 

15 e'i(t,q), au moins pour q compris entre m+1 et m+R, 

(c) et une etape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respectivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'i(t,q), i 

20 etant un indice compris entre 1 et n designant un 
transducteur du reseau d'imagerie. 

2. Procede selon la revendicat ion 1, dans lequel au 
cours de la sous-etape (bll) , lorsqu' au moins certains 
transducteurs (Tl-Tm, T'l-T'm) sont en contact avec un 

25 milieu homogene intermediaire lui-meme en contact avec la 
couche aberratrice, on corrige les reponses impulsionnelles 
hri(t) par repropagation numerique pour simuler des 
transducteurs situes directement au contact de la couche 
aberratrice. 

30 3. Procede selon la revendication 1 ou la 

revendication 2, dans lequel la sous-etape (bl2) comporte 
elle-meme les sous-etapes suivantes : 

(bl21) on determine p matrices de transfert 
H(cok) = [Hri (cok) ] , i allant de 1 a n et r allant de 1 a m, 
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ou Hri«ok) est la valeur, a la frequence cok, de la 
transformee de Fourier de la reponse impulsionnelle hri(t), 

(bl22) on determine pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, n composantes Ei (cok, r) , i variant entre 1 et n, 
telles que F (cok, r) =H (cok) . E (cok, r ) , ou E(cok,r) = [Ei(cok,r)] 
est un vecteur a n composantes, F(cok,r) est un vecteur a m 
composantes Fl(cok,r), 1 variant entre 1 et m, ces m 
composantes Fl(cok,r) correspondant a une focalisation 
souhaitee des ondes acoustiques a la frequence cok sur le 
point de focalisation r correspondant a un transducteur du 
reseau cible, 

(bl23) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, un vecteur de n signaux temporels e (t, r) = [ei ( t, r) ] , 

i variant entre 1 et n, ou ^r^Ei^r).^' en notation 

complexe, ces signaux ei(t,r) etant adaptes pour que 
1' emission de ces derniers respectivement par les 
differents transducteurs i du reseau d' imagerie genere une 
impulsion acoustique focalisee sur le point de focalisation 
r du reseau cible, 

(bl24) une sous-etape de correction des 

aberrations generees par la couche aberratrice entre le 
milieu sensiblement homogene et chaque transducteur cible 
r, ces aberrations etant estimees sur la base des mesures 
precedemment effectuees, les aberrations ainsi estimees 
etant utilisees pour calculer lesdits signaux temporels de 
reference e'i(t,r). 

4. Procede selon la revendicat ion 3, dans lequel au 
'cours de la sous-etape (bl22) on calcule p matrices 
H'Mcok), respectivement par regularisation et inversion des 
matrices de transfert H(cok), et pour chaque transducteur r 
du reseau cible, on calcule le vecteur E(cok,r) par la 
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formule : 

E(tok / r)=H" 1 ((ok) .F(cok, j) . 
.5. Procede selon la revendication 3 ou la 
revendication 4, dans lequel au cours de 1 ' etape <bl22) , 
les composantes Fl(cok,r) du vecteur F(cok,r) correspondant 
a la repartition spatiale du champ desire a la frequence 
(ok, sont egales a 0 pour l*r et egale a 1 pour l=r. 

6. Procede selon 1 • une quelconque des 

revendications 3 a 5, dans lequel au cours de la sous-etape 
(bl24), on assimile la parol aberratrice, au voisinage de 
chaque point de focalisation r correspondant a un 
transducteur du reseau cible, a un f litre a reponse 
impulsionnelle finie, defini a chaque frequence tok par une 
amplitude Gr (cok) et une phase <j> r ((ok), la sous-etape (bl24) 
15 comportant elle-mSme les sous-etapes suivantes : 

(bl241) on calcule, pour chaque frequence 
cok, 1' amplitude Gr(cok) et la phase <|> r (cok), a partir soit 
des signaux ei(t,r), soit des vecteurs E(cok, r) , 

(bl242) on calcule p matrices de transfert 
corrigees H' (cok) = [H' ji (cok) ] , ou 

(bl243) on determine pour chaque 
transducteur r du reseau cible, n composantes E'i(<ok,r), i 
variant entre 1 et n, telles que F (cok, r ) =H' (cok) . E' (cok, r ) , 
ou E'(cok,r)= [E'i(cok,r)] est un vecteur a n composantes, 
F(cok,r) est un vecteur a m composantes Fl(cok,r), 1 variant 
entre 1 et m, ces m composantes Fl(cok,r) correspondant a 
une focalisation souhaitee des ondes acoustiques a la 
frequence cok sur le point de focalisation r correspondant 
a un transducteur du reseau cible, 

(bl244) on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
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cible, un vecteur de n signaux temporels de reference 
e' (t,r)=[e'i(t,r)], i variant entre 1 et n, ou 

e'i(t,r)=£ l E'i((Dk,r). e J °*' en notation complexe. 

7. Procede selon la revendication 6, dans lequel au 
cours de la sous-etape (bl241) , on calcule 1' amplitude 
Gr(cok) et la phase <|> r (cok) comme suit : 



Gr(OJk)=- 



n 

i=i 



V «=i 

^(^0=^|](arg(£'<^,K)))-arg(s(^,r) e - /M ' <)r0,ut )) 



10 ou 



Ei* est la valeur complexe conjuguee de Ei, 
et Ax(rO,r,i)=(d(rO,i)-d(r,i) ) /c, d(r,i) etant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO. 

8. Procede selon la revendication 1 ou la 
revendication 2, dans lequel la sous-etape (bl2) comporte 
elle-meme les sous-etapes suivantes : 

(bl21) on determine p matrices de transfert 
20 H (cok) = [Hri (cok)], i allant de 1 a n et r allant de 1 a m, 
ou Hri«ok) est la valeur, a la frequence cok, de la 
transformee de Fourier de la reponse impulsionnelle hri (t) , 

(bl22') on corrige les matrices de transfert 
H((Dk) pour s'affranchir des aberrations generees par la 
25 parol aberratrice au voisinage de chaque point de 
focalisation r, cette correction etant effectuee a partir 
des reponses impulsionnelles hri (t) precedemment 
determinees, et on obtient ainsi des matrices de transfert 
corrigees H' (oak) , 
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(bl23' ) on determine pour chaque point de 
focalisation r correspondant a un transducteur du reseau 
cible, n composantes E'i(cok,r), i variant entre 1 et n, 
telles que F(cok, r)=H' (ok) . E' <o>k,r) , ou E' <G>k,r) = 
[E'i(cok,r)] est un vecteur a n composantes, F(cok,r) est un 
vecteur a m composantes Fl(©k,r), 1 variant entre 1 et m, 
ces m composantes Fl(cok,r) correspondant a une 
focalisation souhaitee des ondes acoustiques £ la frequence 
(ok sur le point de focalisation r correspondant a un 
transducteur du reseau cible", 

(bl24') on en deduit, pour chaque point de 
focalisation r correspondant a un transducteur du r6seau 
cible, un vecteur de n signaux temporels 

e'(t,r) = [e'i(t r r)] f i variant entre 1 et n, ou 

g, '(^ r )=^^'(^fe,r).g y<ufc# en notation complexe, les signaux 

e'i(t,r) etant lesdits signaux de reference. 

9. Procede selon la revendicat ion 8, dans lequel au 
cours de la sous-etape (bl23' ) on calcule p matrices H'~ 
^cok), respectivement par regularisat ion et inversion des 
matrices de transfert H' (cok) , et pour chaque transducteur 
r du reseau cible, on calcule le vecteur E'(cok,r) par la 
formule : 

E' (oak, r)=H' _1 (cok) . F (cok, j ) . 

10. Proced§ selon la revendication 8 ou la 
revendication 9, dans lequel au cours de l'etape (bl23'), 
les composantes Fl(a>k,r) du vecteur F(a>k,r) correspondant 
a la repartition spatiale du champ desire a la frequence 
©k, sont egales a 0 pour l*r et egale a 1 pour l=r. 

11. Proc6d§ selon I'une quelconque des 
revendications 8 a 10, dans lequel au cours de la sous- 
etape (bl22'), on assimile la paroi aberratrice, au 
voisinage de chaque point de focalisation r correspondant a 
un transducteur du reseau cibie, a un filtre a r6ponse 
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impulsionnelle finie, defini a chaque frequence cok par une 
amplitude Gr((ok) et une phase <t> r (ook), la sous-etape (bl22' ) 
comportant elle-meme les sous-etapes suivantes : 

(bl22'l) on calcule, pour chaque frequence 
cok, l' amplitude Gr(cok) et la phase <t> r (cok) , a partir des 
reponses impulsionnelles precedemment determinees, 

(bl22'2) on calcule p matrices de transfert 
corrigees H' (cok) = [H' ji(cok) ] , ou 

Hj t (a, k ) = H (^)._J_ e -MK) . 

12. Procede selon la revendication 11, dans lequel 
au cours de la sous-etape (bl22'l) on calcule, pour chaque 
frequence cok, 1' amplitude Gr(cok) et la phase <|> r (cok), de la 
maniere suivante : 



<M^=i|:(arg^ ou : 

H*ri designe la valeur complexe conjuguee de Hri f 
et AT(rO,r,i)=(d(rO,i)-d(r,i) ) /c, d(r f i) etant la 
distance entre le transducteur i et le point de 
focalisation r, et d(rO,i) etant la distance entre le 
transducteur i et un point de focalisation particulier rO. 

13. Procede selon 1 1 une quelconque des 

revendications prec§dentes, dans lequel au cours de l'etape 
(c) , on fait suivre la sous-etape (cl) par les sous-etapes 
suivantes : 

(c2) on fait capter par lesdits t ransducteurs du 
reseau d'imagerie des signaux Si(t) ret rodif fuses par le 
milieu heterogene dissipatif, 
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(c3) on convolue le signal de reference emis par 
chaque transducteur du reseau d'imagerie avec le signal 
retrodiffuse capte par ce transducteur, 

(c4) puis on somme les produits de convolution 
ainsi obtenus, 

l'etape (c) dtant renouvelee pour une pluralite de points 
situes dans le milieu sensiblement homogene . 

14. Procede selon I'une quelconque des 
revendications precedentes, dans lequel au cours de la 
sous-6tape (b21), la premiere estimation de chaque signal 
de reference est e'i(t,q) - e' i (ts+0i (q) , qO ) pour chaque 
point - de focalisation q, qO etant I'indice d'un point de 
focalisation proche du point q pour lequel le signal de 
reference a deja ete determine, 9i(q) etant un retard egal 
a une valeur 5i(q)/c, ou c est la celerite moyenne des 
ondes acoustiques dans le milieu, et 8i (q) est egal a une 
difference entre d'une part, une distance entre le 
transducteur i du reseau d' imagerie et le point de 
focalisation qO, et d.' autre part, une distance entre le 
transducteur i du reseau d'imagerie et le point de 
focalisation q, 

15. Procedd selon 1 1 une quelconque des 
revendications precedentes, dans lequel au cours de la 
sous-etape (b2), lorsqu'au moins certains transducteurs 
d'indice v du reseau d'imagerie ne sont pas directement au 
contact de la couche aberratrice, on corrige les signaux 
e' v (t,q) correspondants par repropagat ion numerique pour 
simuler des transducteurs places en contact direct avec la 
couche aberratrice . 

16. Proced<§ selon 1 1 une quelconque des 
revendications precedentes, dans lequel au cours de la 
sous-etape (b22) , on recherche les valeurs a ± (q) et x ± (q) 
pour maximiser le critere de coherence C suivant : 
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C = 





n 


2 


< 


i=\ 


> 



ft 

n -T d <\<x i -g i (t--r i ,q)\ 2 > 



ou 



gi = Si ,(/) ® e i \t,q) , <g> representant l'operation de 
convolution, 

et <> represente une moyenne temporelle . 

17. Procede selon la revendication 16, dans lequel 
au cours de la sous-etape (b22) , les valeurs Xi(q) sont 
calculees en maximisant une fonction d' intercorr61ation, 
pour des transducteurs voisins du reseau d'imagerie, des 
signaux gi(t,q) et g i+1 (t,q). 

18. Procede selon la revendication 16 ou la 
revendication 17, dans lequel au cours de la sous-etape 
(b22), les valeurs ai (q) sont calculees de maniere a 
egaliser sur l'indice i 1' amplitude maximale des fonctions 
5i(t,q) . 

19. Procede selon 1 ' une quelconque des 
revendications precedentes, dans lequel la sous-etape (b22) 
relative a chaque point de focalisation q est realisee 
inunediatement apres la sous-etape (b21) relative au meme 
point de focalisation q. 

20. Procede selon 1 1 une quelconque des 
revendications precedentes, dans lequel le milieu 
heterogene dissipatif est constitue par le cerveau entoure 
par le crane. 

21. Procede selon 1 • une quelconque des 
revendications precedentes, dans lequel : 

soit le reseau d'imagerie et le reseau cible 
sont deux reseaux distincts disposes de part et d' autre du 
milieu heterogene dissipatif, 

- soit tous les transducteurs appartiennent a la 
fois au reseau d'imagerie et au reseau cible. 

22. Procede selon 1 ' une quelconque des 
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revendications precedentes, dans lequel les ondes 
acoustiques sont des ondes ultrasonores . 

23. Dispositif (1) adapte pour mettre en ceuvre un 
precede selon l'une quelconque des revendications 
5 precedentes, ce dispositif comportant un nombre t superieur 
a 2 de transducteurs acoustiques (Tl-Tn, T'l-T'm) destines 
a etre fixes dans des positions predetermines a 
l'exterieur de la couche aberratrice (3), ces transducteurs 
etant commandes par au moins une unite centrale 
10 electronique (CPU) et formant au moins : 

- un reseau d' imagerie (Tl-Tn) qui regroupe un 
nombre n compris entre 1 et t desdits transducteurs, 

- et un reseau cible (T'l-T'm) qui regroupe un 
nombre m compris entre 1 et t desdits transducteurs, 
1' unite centrale electronique etant adaptee pour suivre les 
etapes suivantes : 

(b) une etape d' apprentissage comprenant elle-meme 
les sous etapes suivantes : 

(bl) une sous-etape d' apprentissage de 
20 focalisation du reseau d' imagerie sur le reseau cible, 
etape au cours de laquelle : 

(bll) on determine des reponses 

impulsionnelles hri(t) du milieu heterogene dissipatif, 
respectivement entre chaque transducteur i du reseau 
d' imagerie et plusieurs points de focalisation r situes sur 
la couche aberratrice en correspondance respective avec des 
transducteurs du reseau cible, ces reponses impulsionnelles 
etant memorisees sous forme numerique avec un certain 
echantillonnage temporel qui determine un nombre p de 
composantes f requentielles de la reponse impulsionnelle, de 
frequences respectives cok, i etant un indice compris entre 
1 et n qui designe un transducteur du reseau d' imagerie, r 
etant un indice compris entre 1 et m qui designe un point 
de focalisation correspondant a un transducteur du reseau 
cible et k etant un indice compris entre 1 et p qui designe 
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une composante f requentielle, 

(bl2) a partir de ces reponses 

impulsionnelles, on calcule, pour chaque point de 
focal isat ion r correspondant a un transducteur du reseau 
5 cible, un ensemble de n signaux de reference temporels 
e'i(t,r), i variant entre 1 et n, tels que, si la parol 
aberratrice etait enlevee au voisinage du point de 
focalisation r, 1' emission de ces signaux de reference par 
les differents transducteurs i du reseau d' imagerie 
10 genererait une impulsion acoustique focalisee sur le point 
de focalisation r, 

(b2) une sous-etape de focalisation en un nombre 
R de points de focalisation predetermines situes dans le 
milieu sensiblement homogene, d 1 indices q compris entre m+1 
15 et m+R, cette sous-etape consistant a determiner pour 
chacun de ces points de focalisation q, en s'eloignant pas 
a pas des points de focalisation lam correspondant aux 
transducteurs du reseau cible, des signaux de reference 
e'i(t,q) a faire emettre par les differents transducteurs i 
20 du reseau d' imagerie pour generer un signal predetermine 
focalise sur ledit point de focalisation q, les signaux de 
reference e'i(t,q) etant determines pour chaque point de 
focalisation q en procedant comme suit : 

(b21) une premiere estimation de e'i(t,q), 
25 pour i allant de 1 a 4 f est calculee a partir d'au moins un 
signal de reference e'i(t,q0), q0 etant 1 1 indice d'au moins 
un point de focalisation proche du point de focalisation q 
pour lequel le signal de reference a deja ete determine, ce 
calcul etant fait en utilisant une celerite moyenne des 
30 ondes acoustiques dans le milieu sensiblement homogene (2) , 

(b22) on fait emettre par les transducteurs 
du reseau d' imagerie, par iterations, les estimations 
precedemment obtenues des signaux de reference e'i(t,q), 
puis on capte avec les memes transducteurs des signaux, 
35 Si(t,q) retrodif fuses par le milieu heterogene dissipatif, 
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puis on modifie pour 1' iteration suivante ces signaux de 
reference e'i(t,q) de la maniere suivante : 

ou les valeurs cCi(q) et x ± (q) sont un facteur d' amplitude 
et un retard correctifs, calcules pour maximiser un critere 
de coherence C entre lesdits signaux retrodif fuses, 
lesdites iterations etant arretees lorsque le critere C 
atteint un seuil predetermine, 

(b3) on memorise les signaux de reference 
e'i(t,q), au moins pour q compris entre m+1 et m+R, 

( c .) et un © <§tape de focalisation au cours de 
laquelle, pour au moins un desdits points de focalisation 
q, on fait emettre respectivement par les transducteurs du 
reseau d'imagerie, lesdits signaux de reference e'i(t,q), i 
15 etant un indice compris entre 1 et n designant 
transducteur du reseau d'imagerie. 
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